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Abstract 

The recent positron excess observed in the PAMELA satellite experiment strengthens previous 
experimental findings. We give here an analysis of this excess in the framework of the Stueckelberg 
extension of the standard model which includes an extra U(l)x gauge field and matter in the 
hidden sector. Such matter can produce the right amount of dark matter consistent with the 
WMAP constraints. Assuming the hidden sector matter to be Dirac fermions it is shown that 
their annihilation can produce the positron excess with the right positron energy dependence seen 
in the HEAT, AMS and the PAMELA experiments. The predictions of the p/p flux ratio also fit 
the data. 



1 



Introduction: An excess of the positron flux emanating from the galactic halo has been 
reported by the PAMELA satellite experiment jl| which supports the previous observations 
by the HEAT and AMS experiments 2] but is much more accurate. A remarkable feature 
of the positron spectrum is the turn around and increase in the positron flux with positron 
energy in the range of 10-80 GeV. Here we analyze the possibility that the positron flux 
is arising from the annihilation of particles in the hidden sector. The hidden sector, which 
is defined as a sector with fields which are neutral under the standard model (SM) gauge 
group, has played an increasingly important role since its inception in the formulation of 
supergravity grand unification [3]. Subsequently, it was realized that such hidden sectors 
appear quite naturally in the context of string theory |4|. More recently the hidden sector 
was utilized in the Stueckelberg extension of the standard model. In such an extension, the 
SM gauge group SU(3)c x SU(2)l x U(l)y is supplemented by an extra U(l)x gauge group 
factor where the SM fields are neutral under the U(l)x- However, the mixings between 
the SM and the hidden sector do occur via a connector sector which mixes the gauge fields 
of U(l)y and of U(l)x- Such mixings can occur via the Stueckelberg mass terms[5|] or via 



kinetic mixings 21] or both 




analyzed in 
mixing in 



. Electroweak constraints from the LEP and the Tevatron were 
and an analysis of the dark matter was given with and without kinetic 
lOj consistent with the WMAP 



111 ] constraints. (For related works on the 



and for other works on the hidden sector see 



13] ) . Further recent 



works regarding Stueckelberg extensions in the context of the string and D brane models 
in \\M and for related works see 

n 



151 ]. We note that hidden sectors are also 



can be found 

central to unparticle]16] and ungravity models 12j. 

Positron fraction from annihilation in the hidden sector: The analysis of the positron 
spectrum depends both on the particle physics as well as on the astrophysical models and 
these features have been discussed recently in some detail in [if], Q- Here we focus on the 
fit to the just released data by the PAMELA experiment [l] from annihilation of dark matter 
in the hidden sector in the framework of Stueckelberg extension of the standard model. We 
give now the details of the analysis. 

In general the positron flux arising from the annihilation of dark matter (DM) particles 



is given by 20j, |21 ] 
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where Mp is the mass of the dark matter particle, i] = 1/2(1/4) for the DM particle being 



Majorana or Dirac 



20| | . B is the boost factor which is expected to lie in the range (2-10) 



although significantly larger values have been used in the literature. In the above v e + is the 
positron velocity where v e + ~ c, and p is the local dark matter density in the halo so that p 



lies in the range (0.2 - 0.7)[GeV/cm 3 ] [22j. Further, b(E) in Eq.flU) is given by [23|, [2J, [25] 
b(E) = E (E/E ) 2 /r E , where r E ~ 10 16 [s], with E in [GeV] and E = lGeV. Here (av) H 
is the velocity averaged cross section in the Halo (H) of the galaxy as emphasized by the 
subscript H. In some works (crv) H is replaced by the (av)x } at the freezeout temperature. 
However, such an approximation can lead to significant errors since the ratio (o~v)h/ (<Jv)x f 
can deviate significantly from unity depending on the part of the parameter space one is in. 
The halo function X^^e 1 ) is parametrized as in [21 1. and we consider both the Navarro, Frenk 
and White (NFW) and Moore et. al profiles coupled with various diffusion models. 

The channels that can contribute to the positron flux in the model are DID — > 
ff, W + W~ , . . . where / is any quark or lepton final state of the standard model. In our 
analysis we make use of the dominance of the direct channel poles, specifically of the Z' 
pole, and here the W + W~ final state contribution is much suppressed relative to the // 
final state contribution [6j. We begin by discussing the DD — > ff cross section in Stueck- 
elberg extension of the standard model . Here one has a mixed Higgs mechanism and 
a Stueckelberg mechanism to break the extended SU(2)l x U(1)y x U(1)x electroweak 
symmetry Q]. After such breaking there is the usual mixing between the neutral com- 
ponent of the SU(2)l gauge field A% and the U(l)y gauge field . Thus together the 
Higgs mechanism and the Stueckelberg mechanism generate a 3 x 3 mass matrix for the 
three gauge fields C M , B* 1 , A%. The above model has a non-diagonal mass matrix (M| t ) |5( 
and a non-diagonal kinetic mixing matrix (K.) 7, Q]. A diagonalization of the kinetic term 
can be obtained by a GL(3) transformation. However, such a diagonalization is arbi- 
trary up to an orthogonal transformation. One can choose a combination K' of GL(3) 
and an orthogonal transformation such that Ai 2 = K T Mg t K' has a form similar to that 
of Ml t , i.e., the zero elements of M| t are maintained by the transformation (for details 
see (8|). In the basis where both the kinetic and the mass 2 matrix are diagonal one has 
M. 2 Diag = 1Z T M. 2 1Z. In this basis the interactions between the standard model fermions and 
the vector bosons Z, Z' are given by = f^f [CgP L Z'^ + Cf L P L Z^ + (L -> R)] /, where 
as usual Q = Tl >R +Y LtR /2, T% = 0, P L . R = |(lT7 5 ), and Cf L = T| [g 2 TZ 3i - 7 vT + F7^], 
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FIG. 1: (Left Panel): Positron spectrum including the monochromatic source and continuum flux 
for various halo/diffusion models with (e = 0.006, 5 = 0.00) and p = 0.35 GeV/cm 3 with M z > = 298 



GeV, M D = 150 GeV, t e = 3 x 10 16 s|23j], and B=10; and VLh 2 = 0.13 (calculated by integration 
over the Breit-Wigner pole). Also plotted is the just released PAMELA data |l|], along with the 
AMS-01 and HEAT data [2]. The background flux ratio is the decaying solid (blue) lower curve. 
(Right Panel): An exhibition of the dependence of (crv) on temperature for Stueckelberg models 
as given in the figure with Mp/GeY 6 [150, 153] in steps of 1.5 and Mz> fixed as in the left panel 
of Fig(TTJ). The annihilation near a pole generates a significant enhancement of {ctv)h in the halo 
relative to (crv)x f at freezeout. The natural Breit-Wigner enhancement of {<jv)h obviates the 
necessity of using very large boost factors. 



Cf l R = Q7VI + e 2 "^2i, where Zj = (Z', Z) for i — 1, 2. The photon couples as usual to the 
visible fermionic fields with strength eQ, and 7 is related to e by e = 7(72 /[7 2 + g^} 1 ^- The 
Dirac fermion is assumed to have the interaction gxQxD / -f fJ, DC fl which in the diagonal basis 
gives = DY [C Z D Z>^ + C Z D Z, + (T D Afi D, where CI = gxQxi-ces^ - S s c 0C<j> ), C Z D = 

gxQxiStpC^ + SgS^ty - Ss{s^ - SgCfCy,)), C% = gxQx^C^ - SqS^S^p - Ssic^s^ + SgC^S^p)). 

Here S$ = 5/y/l — 8 2 and the angles 9,(j>,ip appear in the rotation matrix 1Z, and are de- 
fined by tan 9 = r y/g2, tan0 = e = [e — 6][1 — 5 2 }~^ and ip is determined by the relation 
tan 2^ = 2eM ^/M 2 -M^jiM 2 - M 2 + e 2 (M 2 + M 2 - M 2 ,)), where M = vy/g^ + j 2 /2. 
The parameters e and 5 are constrained by the electroweak datajf], Q]. One finds that e and 
5 are both separately constrained so that |e|, \8\ < .06. The action given in [5|] leads to an 
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integrated cross section [101]. 

*// * ££t(i^i 2 + i^i 2 ) • F i + R <atR) ■ f 2 \, (2) 

where F x = 1 + + 4M|,s- 1 (l - 2mJ/s), and F 2 = Sm^ 1 (1 + 2M 2 D /s). Here 

/9/ ijD = (1 — 4m* q/s) 1 / 2 , s = Am 2 D /(l — v 2 /A) and £l,r include the poles 

^ L ' R s s-M 2 z + iT z M z s-M 2 ,+ iT z ,M z , ' 1 ' 

The dominant term in our analysis is the line source arising from the annihilation DD — > 
Z' — > e + e~ and in this case one has ^ F=Final states {crv)p (dN e +/dE') F ~ (av) e + e -6(E' — 
Md) + . . ., where the dots stand for the background terms that contribute to the continuum 
flux. The continuum flux arises mostly from muons and to a much lesser degree from 
taus]28]. Defining Rf as the positron ratio from source / one finds R^/Re is non-negligible 
and decreases with increasing E e + over the dark matter mass (DM) mass range of interest 



and a similar relation holds for the taus[28|. The inclusion of the flux from the continuum 
reduces the needed boost factor slightly, however the line source still dominates at high 
energies. The use of the above in Eq.([T]) yields the primary positron flux $ e + = $ + . We 

(2) 

must add to it the secondary positron flux $ + and then compare it with the electron flux 



($ w + $ w ) . For $^ , ^ , ^ we use the parametrizations of |27l . 1291 ] . For comparison with 
experiment one often defines the positron fraction e + /(e + + e~), and an analysis is given of 
this observable as a function of the positron energy in Fig.([T]) for the Stueckelberg Z' model. 
One finds that the annihilation of Dirac fermions via the Z' pole into e + e~ + fi + n~ gives a 
sufficient kick to generate the necessary turn around in the positron fraction at just about 
the desired value of the positron energy consistent with the relic density constraints. The 
analysis of Fig. (CQ) (left panel) exhibits the theoretical evaluation for several model points. 



Here we consider N 
parametrizations 20 



W min (M2), med and max (Ml) as well as the Moore max (Ml) 



21[. One finds that there is a significant variation in the prediction 



depending on the profile/diffusion model one chooses. However, one finds that the PAMELA 
data does lie in the range of the theoretical predictions. We note in passing that the gamma 
ray spectrum in this model has been discussed in[l(3]. The theoretical predictions cover a 
range which includes the PAMELA data[l|. Further, such a fit determines the dark matter 
fermion mass to be roughly half the Z' mass. 
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FIG. 2: The p/p flux ratio including the TOA correction to the IS spectrum [31[, and with Bp G 
(1 — 10). The green (darker) curves (NFW min) are insensitive to the boost in the ratio, while the 
yellow (lighter) curves (NFW med) allow a boost as large as 5 or even larger. 

In Fig. ([I]) (right panel) we exhibit the dependence of (av) on temperature. The analysis 
shows that (av) can have a significant temperature dependence. Thus the simplifying as- 
sumption often made in assuming that (av) is a constant as one moves from the freezeout 
temperature to the temperature of the galactic halo is erroneous. Specifically the analysis 
shows that the temperature dependence is model dependent and one can generate an en- 
hancement of (av)n in the halo relative to freezeout (av)x f by as much as a factor of 10 
or more depending on the part of the parameter space one is in. Typically the temperature 
dependence is enhanced when the dark matter particles annihilate near a pole from the 
Breit-Wigner which is the case in the analysis here. 

Anti-proton Flux and PAMELA : The p/p flux ratio as recently reported by the 
PAMELA (3^] collaboration indicates a smooth increase with energy up to about 10 GeV 
and then a flattening out in agreement with the background and with previous experiments. 
We note that a suppression of p/p flux ratio is possible in the model presented here. This is 
due in part because the Z' — > W + W~ is suppressed as already discussed. We have carried 
out a detailed analysis of the p/p flux ratio. Our analysis follows closely the work of 2l| 
with fragmentation functions as modeled in Bottino et al and by Bergstrom etal and (p, p) 
backgrounds as in Donato et al and Bringmann et al 31(- The Interstellar (IS) flux has 
been modified for predictions at the Top of the Atmosphere (TOA) which suffers from large 
uncertainties. The results are given in Fig.(T5]) and compared with the recently reported 
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FIG. 3: Fit to the PAMELA and ATIC[3J] data for a heavy Dirac dark matter mass of 688.8 GeV 
with the Breit-Wigner enhancement. The curves in descending order are for the cases for the halo 
profiles listed on the top right hand corner. 

results by the PAMELA collaboration. It is found that the p/p analysis of Fig. ([2]) is fully 
compatible with the recent PAMELA data. It is further observed that the NFW min profile, 
for the p/p predictions, are rather insensitive to a boost factor, while boost factors as large 
as 5 or larger are acceptable in the NFW med model. We note in passing that the p/p flux 
ratio does suffer from larger theoretical uncertainties than the e + /e flux ratio due to a larger 
diffusion length. Further, it is known that local inhomogeneities in the dark matter density 
may lead to very different boost factors for positrons and antiprotons (see, for example, 
Lavalle etal in 3lj|). 

Conclusion: In this work we have shown that the annihilation of the Dirac fermions in 
the hidden sector close to the Z' pole can generate a positron fraction compatible with the 
current PAMELA data. Specifically the model produces the right amount of positron spec- 
trum enhancement with increasing positron energy indicated by the AMS-01 and the HEAT 
data and confirmed by the PAMELA data, and additionally the model can accommodate the 
antiproton constraints. A further support of the model can come from a direct observation 
of the Z' boson at the Large Hadron Collider. 

Note Added: After submission of this paper, the ATIC Collaboration published its data 



on the electron excess 32[]. One can fit both PAMELA and ATIC as well as the p/p flux in 
the model with a change of the dark matter and Z' mass and with a low boost factor of 
10. This fit is given in Fig.Q and includes the continuum flux. Again the Breit-Wigner 
enhancement plays an important role in the analysis. 
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